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Abstract 

We calculate certain string correlation functions, originally introduced as order parame- 
ters in integer spin chains, for the spin-1/2 XXZ Heisenberg chain at zero temperature and 
in the thermodynamic limit. For small distances, we obtain exact results from Bethe Ansatz 
and exact diagonalization, whereas in the large-distance limit, field-theoretical arguments 
yield an asymptotic algebraic decay. We also make contact with two-point spin-correlation 
functions in the asymptotic limit. 

1 Introduction 

Haldane's work [1, 2] on the different ground-state properties of integer-5 and half-integer-S spin 
chains triggered efforts to seek for a quantitative understanding of the gapped ground state of 
integer-S* chains. Among these are the works of den Nijs and Rommelse [3] as well as Oshikawa 
[4], where the following generalized string correlation function was considered: 



0(n,6) = -4(5fexp 



' n-l 
. k=2 



S Z n ■ (LI) 



The authors of [3] introduced Yim. n ^ oo 0{n,n) as an order parameter that characterizes the 
gapped ground state for the 5 = 1 Heisenberg chain and acquires a non-zero value there. 
Kennedy and Tasaki [5] introduced a transformation showing that this is due to a broken hidden 
Z2 x Z2 symmetry of the model. In [4], an attempt was made to generalize this argument to 
integer S > 1 chains. In the same reference, the den Nijs-Rommelse order parameter was 
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considered for 9 / ir. Several subsequent works considered the generalized string correlation 
functions (1.1) for integer spin S > 1 and generic 9, where lim^^oo 0(n, 9) acquires non-zero 
values. The exact calculation for the Valence Bond Solid (VBS) state shows that the correlation 
takes its maximum values near 9 = ir/S [4, 6]. 

Whereas in these works, the focus was mainly on integer spin chains motivated by Haldane's 
conjecture, interest at the same time rose for 0(n,Tr) in half-integer spin chains. Hida [7] 
studied 0(n,ir) for alternating spin-1/2 systems, as this model describes a crossover between 
the gapped 5=1 phase and the isotropic spin-1/2 Heisenberg chain. In that paper, he reported 
the asymptotic form of 0(n,7r) ~ const, n -1 / 4 close to the uniform, isotropic spin-1/2 chain by 
means of a field-theoretical approach (the constant was not known there). This means that the 
string correlation function for the spin-1/2 Heisenberg chain 0(n,ir) decays in an algebraic way 
much slower than the usual spin-spin correlation function. 

Recently, a related string correlation function 



p(n, 9) = ( exp 



k=i 



(1.2) 



was introduced by Lou et al.[8]. They came to the conclusion that asymptotically, for spin 
S = 3/2, 0(n, 9)\ s=3 / 2 ~ — sin 2 (#/2) p(n,9)\ s=1 / 2 - This means that the scaling behaviour of 
p(n, 9)\s=i /2 i s a l so important for S = 3/2, which is supported by the fact that the spin-3/2 
and spin-1/2 chains are considered to belong to the same universality class [9, 10]. Using a 
field-theoretical approach, the authors of [8] found p(n, 9)\ s=1 / 2 ~ const, n 

S 2 /{^ 2 ) W ith an 

unspecified constant, again for the isotropic spin-1/2 chain. 

As far as two-point correlation functions of the spin-1/2 chain are concerned, enormous 
progress has been made in the last decade to obtain exact expressions from Bethe Ansatz 
[11, 12, 13] for short distances [14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28] and from 
field-theory for both the amplitudes and the exponents of the leading terms in the asymptotic 
limit [29, 30, 31, 32]. These results are not restricted to the isotropic point, but cover the critical 
anisotropic regime as well, 



N 

H = JJ2 (SfSf +1 + SfSf +1 + AS?S? +1 ) , (1.3) 
i=i 

with periodic boundary conditions and J > 0. In the following, we use the anisotropy parameter 
7 to parameterize the anisotropy A =: cos 7, with < 7 < ir, such that the isotropic points 
7 = 0, 7r are excluded. 

Given those technical tools from Bethe Ansatz and field theory, in this work we calculate 
p(n,9) and 0(n,9), both for short distances and in the asymptotic limit. We thus obtain the 
exponents and the amplitudes of the leading uniform and alternating parts and verify them by 
the Bethe Ansatz results. Interestingly, the leading asymptotics of the alternating part can be 
directly obtained from those of the uniform part. We furthermore study the limiting values 
9 — > 0, 1 — 7/7r in the asymptotic limit, where contact is made with (SfS^) and (SfS^). 

This article is organized as follows. In the next section, we present the Bethe Ansatz cal- 
culation of p(n,9) and (D(n,9), as well as results from exact diagonalization that we obtained 
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additionally. The third part contains the field-theoretical approach. Numerical comparisons 
between the Bethe Ansatz and field-theoretical results are included in an appendix. Calcula- 
tions not immediately necessary for the understanding of the main text are deferred to further 
appendices. 



2 Exact short distance string correlation functions 

The Hamiltonian (1.3) has been solved exactly by the Bethe Ansatz [11, 12, 13]. In fact the 
eigenfunctions can be constructed in a form of superposition of plane waves, which are called 
the Bethe Ansatz wave functions. The corresponding eigenenergies are obtained by solving the 
Bethe Ansatz equations 

e i fc ,iv = ( _ 1} m-i tt g > +1 2Ae> = 
^ > 11 e i(kj+ki) + i _ 2Ae ifc * 

where M is the number of the down spins. With a solution of the Bethe Ansatz equations (2.4), 
the corresponding eigenenergy is expressed as 

E= JN& +J j2( cos k j -A). (2.5) 

3=1 

Especially the ground state is given by a solution in the sector M = N/2. In the critical region 
— 1 < A = cos 7 < 1, its value per site in the thermodynamic limit N — > oo becomes 



J A Jsin 7 2 f°° dx 
eo = — -— ■ 



f°° dr 

/ (2.6) 

J-oo (cosh 72; — cos 7) cosh irx/2 



Enormous works have been contributed to evaluate the physical quantities of the model based 
on the Bethe Ansatz equations (2.4) [13]. They, however, are usually limited to the bulk 
quantities. Especially, the exact calculation of correlation functions still is a difficult prob- 
lem. Only for A = 0, where the system reduces to a lattice free-fermion model after a Jordan- 
Wigner transformation, arbitrary correlation functions can be calculated by means of Wick's 
theorem [35, 36]. Especially, the two-point spin-spin correlation function is simply given by 

(5;s; + ,) = -(i-(-if)/(2A 2 ). 

There have been many attempts to evaluate the correlation functions for general A. However, 
explicit exact evaluations of the correlation functions have become attainable only recently. 
For example, the following exact values for the spin-spin correlation functions (SjSj +k ), were 
obtained up to k = 7 for A = 1 [24] and up to k = 8 for A = 1/2 [37]: 
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• A = 1 



(SjSj +2 ) 
(SjSj +4 ) 



(S*S Z 



j i+6 
(SjSj +7 



• A = 1/2 



(SjSj +1 

I QZ QZ 
I QZ QZ 
I QZ QZ 

(SjSj +5 
(SjSj +6 

I QZ QZ 

(SjSj +8 



— - - In 2 = -0.147715726853315 
12 3 

— - - In 2 + -C(3) = 0.060679769956435 
12 3 4 SW 

1 37 14 3 

— - 3 In 2 + -C(3) - y In 2 • C(3) - -C(3) 2 

125 25 
-— C(5) + y In 2 • C(5) = -0.0502486272572352 • 

i L-H ln2 + ^ C(3 )_541»2. < (3,-^ 



875 w . 1450, . . 275 „, , , N 1875 , 



., In2-C(5)-— C(3)-C(5)-^C(5) i 



3185 >/w , 1715, 6615 „, . 

f — C(7) - -r- In 2 ■ C(7) + -^-C(3) • C(7) 



64 ' w 4 
0.0346527769827281 • • • 

-0.0308903666476093- 

0.0244467383279589 • • • 

-0.0224982227633722 • 



-- = -0.125, 
7 



32 



(2.7) 



256 



0.02734375, 



401 



= 0.02447509765625, 



0.0109942555427551 



-0.0110844789305701' 



16384 
184453 

16777216 
95214949 

8589934592 

1758750082939 

— — — = 0.0062483354772489 

281474976710656 

30283610739677093 

= -0.0065667113109326' 

4611686018427387904 

5020218849740515343761 
1208925819614629174706176 



0.0041526277032786- 



(2.8) 



Here C(2& + 1) is the Riemann's zeta function at odd arguments. Note that the nearest- 
neighbour correlation function (SjSj +1 ) can be derived immediately from the ground state 
energy (2.6). So these values have been known long before. We also remark (5|5| +2 ) for A = 1 
was obtained some decades ago by Takahashi [38] by his ingenious study of the half-filled Hub- 
bard chain in the strong coupling limit. Other results are due to recent developments of the 
study of the correlation functions. Note that even for general A, the exact analytic expressions 
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have been obtained up to k = 3 [21]. Such progress has enabled comparison with the field- 
theoretical prediction of the asymptotic behaviour as well as other numerical methods such as 
numerical diagonalization. 

It is interesting to note that the calculation of the spin-spin correlation functions Eq. (2.7) 
and Eq. (2.8) rely on the generating function, defined by 



p .-(n{0+ s i)+"(i -*?)})■ (2 - 9) 



Here k is a parameter. Once the generating function Eq. (2.9) is calculated, the two-point 
spin-spin correlation function can be obtained by the formula 



i d 2 



i 



{SfS ^ = 2d^ ^ " + P ^ K =i " I" (2 ' 10) 

The generating function (2.9) together with its relation to the two-point spin-spin correlation 
function (2.10) was introduced by Izergin and Korepin [39, 40] (see also the book [12]). Subse- 
quently it was utilized to discuss a certain long-distance asymptotic behaviour [41, 42] as well as 
to obtain several different forms of multiple integral formulas [17, 18]. However, it was only quite 
recently the generating function P£ was explicitly calculated for A / 0, namely, up to n = 8 for 
A = 1 [24] and up to n = 9 for A = 1/2 [37]. 

Now one will readily find P£, Eq. (2.9) and the string correlation function p(n,6), Eq. (1.2) 
are connected as 

p(n,9) = K'^P: (2.11) 
K=e~ ,e 

Then we can calculate some exact values of p(n,6) for A = 1 and A = 1/2. Moreover, since the 
generalized string correlation function O(n,0) (1.1) is related to p(n,9) as 



0{n,9) 



sin 2 f 



p(n, 9) - 2 ( cos - J p(n - 1, 9) + ( cos 2 ) p{n - 2, 9) 



(2.12) 



we can also evaluate the generalized correlation functions for A = 1 and A = 1/2 (cf. Appendix 
A). They are plotted in Fig. 2.1 and Fig. 2.2. From the figures one observes: 

• For even n (> 4), 0(n,9) is always positive with a period 2ir. It has a single max- 
imum at 9 = 7r and a minimum at 9 = 0. Recall that 0(n,ir) = (2i) n (11^=1 ^k) anc ^ 
O(n,0) = -4(Sf^). 

• For odd n, 0(n,9) has a rather complicated structure with a period 4ir. In this case, 
0(n, 7r) and C(n,37r) are always zero as they should be. 

We give some exact values of 0(n, it) = (2i)™(nfe=i <S|) for A = 1 and A = 1/2 in the following: 
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Figure 2.1: 0(n,0) for A = 1 
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Figure 2.2: 0(n,6) for A = 1/2 
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• A = 1 



1 4 

0(2,71") = — + -ln2 = 0.5908629074132604- •• 
3 3 

0(4, t) = ^-y ln2 + ^(3)-f ln2- C (3)-|c(3) 2 

-— C(5) + — In 2 • C(5) = 0.4914453923615522 • • • 
6 3 

0(6,71-) = 0.4403016697026268 •• • 

0(8,71") = 0.4072424147596208 •• • (2.13) 



A = 1/2 



0(2, 7T) 



1 

2 

1595 



0(4, tt) = = 0.389404296875 

v ' ; 4096 

719423395 

0(6, tt) = — — — = 0.3350076242350041 • • • 

v ' ' 2147483648 

346891287109196331 

0(8, tt = = 0.3008802296800668 •• • (2.14) 

V ' ' 1152921504606846976 V ; 

One observes that 0(n,7r) for n=even decays very slowly as n increases. Namely as mentioned 
in the introduction, for A = 1, the asymptotic decay 0(n, tt) ~ n -1 / 4 was given by Hida [7] and 
more generally 

0(n,<9) ~ n ^ (2.15) 

by Lou [8]. In the next section, we shall both generalize this asymptotic formula to the more 
general — 1 < A < 1 case and determine the amplitude by making use of field theory. Further- 
more since the formula Eq. (2.15) does not explain the difference of the periodicity with respect 
to the parity of n, we shall consider some subleading terms more carefully. We remark that 
p(n,0), Eq. (1.2), shares periodicity properties analogous to O(n,0). In fact it is easy to see 
that p(n, 0) is expanded as 

n 

p(n,8) = £p nJ cos [g -j) 0\ , (2.16) 

where the coefficients P n j are the summation of the diagonal density matrix elements in the 
sector with j down spins. Note that P n j = P n>n -j. From Eq.(2.16), one can immediately find 

p(n,0 + 2*) = (-l) n p(n,0) (2.17) 

• A = case : 
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Let us comment on the string correlation functions for A = 0. In this case, a simple deter- 
minant formula for the string correlation function p(n,0) exists (c.f. Ref.[43]) 



p(n, 6) = det 



cos- 8 jtk + ( isin- ) M j)k 



M. 



j,k=i 



-1)" 



j-fc+i 



, 7T j — k 



j — k : even 
j — k : odd 



(2.18) 



Using this formula one can evaluate the exact numerical values up to the order of n ~ 10000 
easily, for example, by Mathematica on a standard PC. We give the exact values in Table 2.1 up 
to n = 1000. 





n 


5 


10 


20 


50 


100 


200 


500 


1000 


p(n 
p(n 
p(n, 


.tt/4) 
.tt/2) 
3tt/4) 


0.884857 
0.605461 
0.289075 


0.866761 
0.564975 
0.291125 


0.848076 
0.516684 
0.234367 


0.824090 
0.459997 
0.177503 


0.806421 
0.421580 
0.144555 


0.789137 
0.386481 
0.118072 


0.766860 
0.344597 
0.0906455 


0.750428 
0.315979 
0.0743400 



Table 2.1: Exact values of p(n,0) for A = 
This determinant is also expressed as a Toeplitz determinant 



<V if? 



p(n,6) 



e^det 



27r A) J j,fc=i 



1 



-W 



71 

:0<q<- 

ir 3ir 
'■ tt < 0. < "tt 

& I 
: — < q < 2tt 
2 



(2.19) 

There are some mathematical results known on the asymptotic behaviours of Toeplitz deter- 
minants as n — ► oo. Assume 9 / 0,27r, then "the generating function" a(q) of the Toeplitz 
determinant has jump singularities at q = tt/2 and q = 3ir/2. In such a case, we can invoke the 
(generalized) Fisher-Hartwig conjecture [44, 45], which brings about an asymptotic formula for 
< 6 < 2n as 



p(n,0) 



P { AL(n,9) + (-irpW ym (n,e), 



n 



-2(-A +fe) 2 4 _ ( _|_ +fe)2 



1 2 



0(1+ * )0( i_ +*) 



(2.20) 



Here G(z) is the Barnes G-function defined by 



<?(,+ !)- exp (-1, - \ b + 1)/) { (l + |)'exp + } 

where 7 = 0.5772156649 •• • is the Euler-Mascheroni constant. Each term PAsym( n '^) decays 



algebraically with the exponent — 2 (— ^ + ky . Then the dominant term is PAs ym ( n '^) f° r 
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n 


5 


10 


20 


50 


100 


200 


500 


1000 


PAsym(n,ir/4) 

PAsym(«, 7T/2) 
/9Asym(n, 37r/4) 


0.884970 
0.605720 
0.289328 


0.866783 
0.565076 
0.291316 


0.848081 
0.516705 
0.234403 


0.824091 
0.460000 
0.177507 


0.806421 
0.421581 
0.144555 


0.789137 
0.386481 
0.118072 


0.766860 
0.344597 
0.0906455 


0.750428 
0.315979 
0.0743400 



Table 2.2: Numerical values obtained from the asymptotic formula pAsym("-, 0) for A = 

< 6 < 7T, and (-l) n PAsym( n ' 0) for 7r < < 27T. We refer the reader also to Refs. [46, 47] for 
more information about the (generalized) Fisher-Hartwig conjectures. 

Numerical values calculated from Eq. (2.20) are listed in Table 2.2. Good agreement is found 
with the data in Table 2.1. In this context, it is remarkable that they coincide within at least 
three digits even for small sizes as n = 10. Finally let us note a further exact result for p(n, ir) 
at A = 0. Since p(2m — l,7r) = 0, we consider p(2m,ir), which is given more explicitly as 



' 2m 



p(2m,n) 



(-l)^ ^11 SjJ = ("l) m det [M hk Y™ =l 

2m m m , ' — k \ 2 m 



exp 



i r 

.~2i 



dt 1 



-mt 



com 1 I 2 exp 



t cosh 2 (t/ 4) 

'°° dt 



i r c 
2/0 



e 'tanh 2 , 
2 L t \ Am 



com 



-1/2 



1 m H m 

32 2048 



379 
65536 



m" u + 



(2-21) 



where 



co 



exp 



1 



dt 



e" 4 '- 



1 



cosh 2 1 



Oil 



.2 Jo * 

Here we have used an integral formula for the logarithm of the Euler gamma function, 



io g r 



(*)= / 
j 



(z-l) 



1 _ e (-^+l)t 



—dt, (»(z) > 0). 



(2.22) 



Thus we can obtain the asymptotic expansion to an arbitrary order in this case. Note that the 
leading term is consistent with the formula Eq. (2.20) with 9 = it. 



3 Asymptotic behaviour of string correlation functions 

In this section, we will discuss the asymptotic behaviour of the string correlation functions for 
the critical region — 1 < A < 1 (that is ir > 7 > 0) by use of field theoretical arguments. Thus 
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the aim is to find coefficients Dj and exponents uj such that 
p{n,9)-Y^Dj{9 n )n-"^) 



lim 

n— >oo 



: An (0,7) (finite), m = l,2,.... 



(3.1) 



n -f m (0,7) 

The exponents are increasing with j, i.e. i/j < Vj+i- The amplitudes and exponents depend 
on the parameters 0,7 of the model and of the function p. Instead of Eq. (3.1), we use the 
shorthand notation 



p(n,0)~^^(0, 7 )n-^ 



(<?.7) 



The important point to remember is that the asymptotic expansion is defined in the limit n — > oo. 

We first present the results obtained so far within the field-theoretical framework and give 
details of the derivation in the following section. 



3.1 Results 

We find the following asymptotic expansion of the string correlation function for < < it: 



p(n, 0) = ( exp 



k=l 



D(0, j)n- Ul ^ (l + O (n"*W)) 
+(-l) n L>(27r - 0, 7 )n^ l( "' 7) (l + O 



n 



-5( 7 ) 



(3.2) 



with 



"1(0,-7) 
5( 7 ) 

D(6, 7 ) 



e 2 vr 



47T 2 7T — 7 



7T 



-4 



7T — 7 



r 



2-2r] 



2-2-q 
oo 



20Fr' 

x exp 



e 2 /(4??7r 2 ) 



sinh 



2tt 



? 2 e" 2 * \ dt 



r 



\ sinht cosh(l — rj)t sinhr/i 4rjn 2 J t 

\1-2ttK 2 ) 



(3.3) 



x exp 



2-4ttH 2 
oo / 



sinh 2 ^-t 

2,7V 







2 e~ 2 * 



dt 

sinht cosh(l - 2vri? 2 )t sinh27ri? 2 i V 27r -R J 2vr / T 



(3.4) 
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In Eq. (3.3), Lukyanov's notation is used with 77 



7T — 7 



whereas in Eq. (3.4), the anisotropy is 



written in terms of the compactification radius R with 2ttR = r\. 

Since p(n,9) = p(n,—9), the result (3.2) is readily extended to the domain —tt < 9 < 0. 
Thus p(n,9) is known in the fundamental domain — ir < 9 < it (note p(n, 9 = 0) = 1, trivially). 
The periodicity Eq. (2.17) then yields p for all values of 9. 

We note the following limiting values of the coefficient D(9,j): 

• D(9 = 2tt7],^) = 2(1 — r}) 2 A, where A is the coefficient of the leading term in an asymp- 
totic expansion of the uniform part of (crfa*), namely: (af [32]. Then, 
as vi(9 = 2ttt], 7) = rj, we have the asymptotic equality (note that 1 — rj = 7/71") for 

7r/2 < 7 < 7T 



p(n,2^)~2(^) 2 (af<) u , 7r/2<7<vr 



(3.5) 



for the leading order of the uniform part (in order to facilitate comparison with Lukyanov's 
results, we use the Pauli-matrices a u = 2S^). This is checked immediately at the free 
fermion point A = from Eq. (2.21). 



-D(0,7) = 1, whereas 
Jim 0(2* -0,7)^ 



A* 
2 



7T^ 



2^r 



l/r, 



x exp 



2-2)? J _ 

sinh((2r/ - l)t) 



U0 



27/ - 1 



sinh^t) cosh((l — rj)t) rj 



e~ 2t I- 



dt 
1 



(3.6) 



The last equation is proved in Appendix B. Following Lukyanov's notation [32], A z denotes 
the coefficient of the leading contribution in the alternating part (afa^) a of the a z -a z - 
corr elation function, namely 



(°"l°"n)a 



(-1) 



n— 1 



n 



i/ v 



(3.7) 



In order to obtain the asymptotics of the generalized string correlation function, we first express 
it in terms of p(n,9) according to Eq. (2.12). Then, using the above results, the asymptotic 
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behavior of 0(n, 0) is obtained for < < it: 



O(n,0) = -4( Sfexp 



' n-l 

k=2 



QZ 



D{0^)n- v ^ 6 ^ 



2 cosf i^(fl, 7 ) cos f(2 cos 1-1)^(^,7) (^(^7) + !) 
Tan — z~pr 1 7r-z k r 



+ (-l) n D(2ir -0,j)n 



4 cos 2 | n 

-J/l(27T-fl,7) 



sin 2 | 



2 (9 cosf i^(27r- (9,7) cos |(2cos § + 1) Vl {2ic - 0,7)(i/i(27r - 0,j) + 1) 



cor - + 



4 sin 



2 
4 



n 



+ 



sin 



+ 



(3. 



Let us consider the limit — ► of the asymptotic formula Eq.(3.8). The first two terms of the 
uniform part in (3.8) vanish in this limit and the third term gives 



limD(0,-f)n- Ul( - 9 ^ 
9^0 



cos§(2cos§-l)^(0,7)M0,7) + l) 



sin 



2 e 



n- 



ir(ir — 7)n 2 Tx 2 r\n 2 



(3.9) 



Because of the relation Eq. (3.6), the leading alternating part yields 

hm(-irZ)(27r-0,7)n-^( M ^cot 2 ? = tllAl 



(3.10) 



In order to get the correct leading alternating term in the limit — ► 0, we should also consider 
the leading alternating part for — it < < 0, which reads 



(-l) n D(27T + 0, 7 ) n -"i(2T+o,7) 



^ 2 fl < COS | |^(27T + g,7) , 
cot t" /) t 



sin 



n 



(3.11) 



in addition to (3.8). Then we have similarly to Eq. (3.10) 



lim(-l) n I>(27r + e,7)n-^ (2,r+fl ' 7) cot 2 °- = 



(3.12) 



in the limit 0^0. Collecting the terms (3.9), (3.10) and (3.12) yields 

limO(n,fl) ~ 1 + 



0->o " ' ' n 1 /*? ' TT 2 7]n 2 ' 



(3.13) 



Eq. (3.13) coincides with the leading asymptotic behaviour of two-point correlation function 
-4 (StS% =-(*&). 
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3.2 Derivation 



3.2.1 The string correlation function p(n,0) 

An effective field theory describing the low-energy-excitations of the XXZ-chain in the critical 
regime < 7 < ir has been derived by Lukyanov [31]. At zero magnetic field, the corresponding 
Hamiltonian density H reads 

+e 2 [\ + J 2 L (x)J 2 R (x) + A_ [Ji(x) + J%(x)]] +..., (3.14) 

where the dots denote terms with scaling dimensions higher than those given explicitly. The 
dimensionless constants v, A, A± are known exactly from Bethe Ansatz [31]. The lattice constant 
e has the dimension of a length, whereas the dimension of TL is 1/length 2 . Above, TL is the 
sum of a Gaussian model and irrelevant operators, the latter with scaling dimensions \/(ttR 2 ) 
and 4. The Hamiltonian is written in terms of the bosonic field cp and its momentum II, where 
[ip(x) ,H(y)] = iS(x — y). The left- and right current operators are then defined as Jlr{x) = 
(U(x) ± <p{x)). 

Within the same approach, the effective S z -operator reads 

S* ee S'(x) ~ ^9Mx) + E (-l) m e^ 2 ^C m cos (^ 1 V>(x 

+descendants , (3.15) 

where x = ej. The constants C m have been determined in [32]. We do not write down the 
descendant fields here explicitly, but only note that if a primary field has scaling dimension A, 
then the descendant fields have scaling dimension A + £, where I is a certain positive integer. 

To arrive at an asymptotic expression for p(n, 9), we first apply the Euler-MacLaurin formula 
to the sum in the exponent: 

X> 2 = e- 1 [ S*{x')Ax'- l -[S z {Q>)-S z {x)] + 0{d»S z ) 
k=i Jo 

= ^ (P(x) - <p{0)) + 0( £ ( 2 ™ +1 ) 2 /(^ 2 )+m £ »), (3.16) 

(p > 1 integer) from which one concludes that the only cutoff-independent contribution in the 
integral stems from the first term in the last equation. We expand the exponent with respect to 
e and arrive at 

p(n, d) ~ (exp i-^{<p{x) ~ <P(0)) (l + 0(£ 2k ^ 2rn+l ^ , e 2fc ^)) ^ (3.17) 

where the positive integer k originates in the expansion of the exponential function. From this 
we conclude that the leading exponent of the uniform part is z/i(6>, 7), and the leading Euler- 
MacLaurin corrections to this have exponents ^(#,7) + 2, ^1(0,7) + l/(27ri? 2 ). 
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Thus in order to determine the amplitude of the leading term, we have to calculate 
(exp [^2^r( 1 p( x ) — 9?(0))] ). In the field-theoretical setting of massless Bose fields considered 
here, this quantity is defined only up to a multiplicative constant A with dimension 1 /length 
[32]. It has become custom to choose it such that ("CFT normalization condition") 

A " 2 /(2-)( exp [ia^x) - 99(0))]) = \x\- a2 / (2w \ (3.18) 

This means that we have to introduce a constant D(9,j) as follows: 

£>(0, 7 ) 



p(n,0)~ L>(0, 7 )(exp 



2ttR 



(<p(x) - <p(0)) 



n (e/2nR)y(2n) 



(3.19) 



for the leading decay of the uniform part. Because of the symmetry p(n,—Q) = p(n,9), this 
result is valid for — it < 6 < it. Let us defer the calculation of the coefficient D to the next 
paragraph and first determine the leading exponent of the alternating part. This is obtained 
directly by exploiting the periodicity Eq. (2.17). Together with Eq. (3.19), it implies that 

pin, ir<9<3ir) ~ D(9 - 2ir, 7 )„-^(*- 2w .t) (3.20) 
p(n, -3ir<9< -it) ~ D(9 + 2ir, 7)7^ ( e+2 ^). (3.21) 

The exponents are expected to depend continuously on the parameters 9, 7. Thus for < 9 < it 
(—-it < 9 < 0), Eq. (3.20) (Eq. (3.21)) yields the leading contribution to the alternating part, 
which is next-leading with respect to the leading decay given in Eq. (3.19). 

What are the exponents of the next-leading contributions? In Eq. (3.19) we have tacitly 
assumed that the expectation value is taken with respect to the unperturbed Gaussian part 
of the Hamiltonian (3.14). However, there are additional contributions in Eq. (3.14), with 
scaling dimensions A = l/(7ri? 2 ),4. As argued in [34], they lead to exponents ^1(0,7) + k(A — 
2) in (exp [^2^n( ( p( x ) ~ ^(O))]), where the integer k denotes the order of the perturbational 
expansion. Since there is no contribution of the cos-operator for k = 1, the next-leading exponent 
stemming from this contribution is 1^(0, 7) = ^1(0,7) + £(7) with 5(7) = 4tt/(tt — 7) — 4. On 
the other hand, the first-order contribution of the A-t-operators yields an exponent fi(9, 7) + 2. 
This latter one is always larger than 1/1(0,7), ^i(0 — 2-zr, 7) (for < 9 < ir) and we discard it here. 
Thus 1^2(9, 7) yields the next-leading exponent in the uniform part. According to the periodicity 
argument, the next-leading exponent in the alternating part for < 9 < it is then v<i(9 — 2tt, 7). 

We now focus on the coefficient D(9,j). The result given in Eq. (3.3) is a conjecture based 
on the work [33]. The following tests of this conjecture have been performed: 

• For 7 = it/2, one can show that D(9,tt/2) reduces to (2.20), namely 



D(9,tt/2) = 4~3^ 

This equality can be checked by means of an integral representation of the Barnes G- 
function (see Appendix C). 



G\ 1 + £-)G[ 1- 

2tt J \ 2vr 



(3.22) 
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• Numerical comparisons for A = 1/2 between the exact data from the Bethe Ansatz (for n = 
9) and the asymptotic results (3. 2), (3. 8) have been performed for 9 = 7r/4, tt/2, 37r/4, it. 
In all cases, very good agreement is found. Similarly, we compared with the data obtained 
by the numerical diagonalization up to a system size of N = 28 lattice sites for general A 
(see Appendix A). 

Our conjecture for D is based on arguments similar to the conjecture for the coefficient of the 
leading decay of (a'fa^), cf. [32], [33]. In [33], the expectation value of (exp [im?]) in a massive 
Sine-Gordon model with an operator cos(/3i9) is determined, 

(exp [icn?]) = (£m) fl2/(4T) N(a, 0), (3.23) 

where m is the particle mass associated with the field Since in that problem, a x ~ e l2,r 
with an a priori unknown amplitude, calculating the amplitude of the leading decay of (cf a^) 
with respect to a Sine-Gordon model with an operator cos(/3$) is very similar to our problem of 
determining D. 

An explicit value for N(a, (3) in Eq. (3.23) is conjectured and confirmed explicitly in certain 
limiting cases in [33]. The authors then calculate (afcr*) ~ y ^( r ?)-^(l/ r ?) 2/R)n~ ri by making use 
of the fact that this correlation function is known explicitly for the massive XYZ-model close to 
the critical XXZ-point, namely (afafj) m ~ A m (em)~' n with a known coefficient A m . This allows 
for the deduction of A{rj). 

In our case, the field tp is related to i? by dtp = d x $. However, the problem of calculating 
D is completely analogous to the calculation of A{rf) sketched above, with a Sine-Gordon-term 
cos(2ip/R) in the Hamiltonian. The only unknown is the string function p m {6) in the massive 
XYZ-regime. We know that 

p m (9) = C m (R)(emf/^ 2 /^ (3.24) 

with an unknown constant C rn depending on R. Note that in the massive regime, we cannot 
relate 2/R to 7, but rather take it as the constant in the Sine-Gordon-term cos(2p/R). On the 
other hand, the results in [33] tell us that 

p m (9) = D N(9/(2ttR), 2/R){em) {e /^ R ^ 2 l^ (3.25) 

with a known coefficient N(9/(2ttR), 2/R). By comparing Eq. (3.25) with Eq. (3.24), one obtains 
D in terms of 9, R and the unknown C m (R). We find that C m (R) = 2(1 - r?) 2 = 2(1 - 2ttR 2 ) 2 
yields excellent agreement with the numerical data as described above. This results in the 
coefficient D(9,j) as given in Eqs. (3.3), (3.4). 

We finally comment on the isotropic case, 7 = 0. Here, 1^1(^,7 = 0) = # 2 /(47r 2 ), in agreement 
with the result of [8]. However, we expect that a logarithmic dependence of the amplitude on 
the distance occurs, similarly to what happens for the two-point functions [30, 31, 32]. We leave 
the study of this case as a project for future research. 

3.2.2 The generalized string correlation function 0(n,9) 

From Eq. (2.12), the asymptotics of 0(n, 9) is obtained once the asymptotics of p(n, 9) is known. 
It is nevertheless instructive to perform a consistency check of this result by calculating the 
asymptotics of O directly by using field-theoretical arguments. 
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Therefore, one might be tempted to take the asymptotic expansion of S z (x), Eq. (3.15), and 
insert it into Eq. (1.1). However, in such a calculation the leading terms given in Eq. (3.2) would 
be absent. We are thus lead to use the following asymptotic expansion for the S^-operators at 
sites 1 and n involved in 0(n, / 0): 

S z (x) ~ s + ^j^^ix) 

+ jr (-l)%( 2 ™+D 2 /(4^ 2 ) Cm cos f^+l^ x )\ + descendants , (3.26) 

m=0 ^ ' 

with x = ej. The asymptotic expansion starts with a finite constant so- For the asymptotics 
of the phase factor in O(n,0), we still use Eq. (3.15). Carrying out the same calculations as 
above, one finds Sq = tan 2 (0/4), which vanishes for 8 = 0. The intriguing point is that we have 
to modify the asymptotic expansion for the spins at sites 1 and n in O(n,0) without modifying 
the Hamiltonian. Namely, it looks as if in the asymptotic limit, the phase operator in 0(n,8) 
acts as a local field on the edge spins. 

4 Conclusion and outlook 

We evaluated the string correlation functions p(n,9) and O(n,0) for the critical anisotropic 
spin-1/2 chain. For small n, exact results were obtained from the Bethe Ansatz, whereas in 
the asymptotic limit, both the amplitudes and the exponents of the leading decay could be 
determined from field theory. The field-theoretical results agree well with the Bethe Ansatz 
data. Especially, for A = 0, the asymptotics could be confirmed directly from the Bethe Ansatz 
results. 

Most interestingly, the leading decay of the two-point xx-correlation function was recovered, 
Eq. (3.5). Whether this result has a physical background has to be clarified. As far as the limit 
— ► in O(n,0) is concerned, we have recovered the expected result Eq. (3.13). However, the 
rather heuristic expansion (3.26) in the field-theory for 6 / deserves further investigations in 
the future. 
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Appendix A Numerical values of string correlation functions 

For A = 1 and 1/2, the string correlation functions p(n,6) and O(n,0) can be evaluated ana- 
lytically up to n = 8 and n = 9, respectively. Here firstly, we list their precise numerical values 
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for 9 = 7r/4, 7r/2, 37r/4, 7r up to n = 8, based on these analytical expressions (see Table A.l and 
Table A. 2). Note that p(l,6) = cos | irrespective of A and therefore we have 

p(l,|) = cos I = 0.923880, p (l, |) = cos j = 0.707107, 

/ 37T \ 37T 

p i— = cos — = 0.382683, p(l,vr) = 0. 



4 / 

Note also that 0(2,6*) = — 4(SfSf) irrespective of 6> by its definition. 



n 


2 


3 


4 


5 


6 


7 


8 


p(n,7r/4) 


0.940083 


0.915627 


0.925111 


0.910171 


0.917092 


0.90616 


0.911707 


p(n,n/2) 


0.795431 


0.685542 


0.744898 


0.671293 


0.718266 


0.66085 


0.70065 


p(n, 3tt/4) 


0.65078 


0.362761 


0.565509 


0.349604 


0.521325 


0.339972 


0.492564 


p(n,7r) 


0.590863 





0.491445 





0.440302 





0.407242 


C(n,7r/4) 


0.590863 


-0.224243 


0.24353 


-0.120692 


0.170343 


-0.0786391 


0.137344 


0(n,7r/2) 


0.590863 


-0.171628 


0.34622 


-0.0787594 


0.282725 


-0.038569 


0.2504 


C(n,37r/4) 


0.590863 


-0.0928846 


0.44891 


-0.0352563 


0.394313 


-0.00917944 


0.361674 



Table A.l: Exact values of string correlation functions for A = 1 



n2345 6 7 8 



p(n,7r/4) 


0.926777 


0.909081 


0.909299 


0.900034 


0.899811 


0.893662 


0.893337 


p(n, it/2) 


0.75 


0.668437 


0.692139 


0.644847 


0.661899 


0.628299 


0.641883 


p(n,3w/A) 


0.573223 


0.346957 


0.477542 


0.32521 


0.429591 


0.310018 


0.398987 


p(n,ir) 


0.5 


0. 


0.389404 


0. 


0.335008 


0. 


0.30088 


0{n, tt/4) 


0.5 


-0.101049 


0.14059 


-0.0285505 


0.0880858 


-0.00509166 


0.0689805 


0(n, it/2) 


0.5 


-0.0773398 


0.243652 


0.000467493 


0.192033 


0.0293118 


0.168568 


C(n,37r/4) 


0.5 


-0.041856 


0.346715 


0.012332 


0.29362 


0.033798 


0.263161 



Table A. 2: Exact values of string correlation functions for A = 1/2 



For A = 1/2 let us compare the results above with the numerical value of the asymptotic 
formulae Eq. (3.2) and Eq. (3.8) with 7 = vr/3 in Table A.3. 



n 


2 


3 


4 


5 


6 


7 


8 


PAsym(™,7T/4) 


0.926694 


0.909865 


0.909106 


0.900388 


0.899692 


0.893859 


0.893259 


PAsym(n,7T/ 2 ) 


0.751733 


0.669839 


0.692208 


0.645454 


0.661862 


0.628622 


0.641843 


PAsym(n, 37r/4) 


0.577912 


0.348016 


0.478151 


0.325667 


0.429729 


0.310258 


0.399020 


PAsym(n,TT) 


0.506119 





0.390271 





0.335222 





0.300940 


As ym(n,7T/4) 


0.306262 


-0.0812361 


0.116753 


-0.0204905 


0.0793274 


-0.000774240 


0.0644446 


As ym(n,7T/2) 


0.402260 


-0.0722680 


0.233892 


0.00275605 


0.1888360 


0.0305388 


0.167037 


e>Asym(n,37T/4) 


0.477679 


-0.0413030 


0.344917 


0.0128694 


0.293024 


0.0341302 


0.262857 



Table A. 3: Asymptotic formulas of string correlation functions for A = 1/2 
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We find the exact values and the asymptotic formulas are in good agreement especially for 
p(n,6). The deviation is somewhat larger for 0(n,9), for which we probably need higher order 
corrections to the asymptotic formulas to achieve better agreement. 

To confirm our asymptotic formula further, we have calculated p(n, 9) numerically for sev- 
eral values of A(= ±0.3, ±0.7, —0.5) by means of the exact diagonalization for finite systems 
N = 20 ~ 28. Then we have applied an extrapolation according to cq + c\/N 2 + C2/N^ + C3/N 4 + 
C4/N 5 and estimated pNum(n,^) in the thermodynamic limit. These values are compared with 
our asymptotic formula pAsym(n, 9) in Tables A.4-A.8. We conclude that our asymptotic formula 
gives fairly precise values for all ranges of A in the critical region. 



n 


2 


3 


4 


5 


6 


7 


8 


PNum(«-,7r/4) 

PAsym(n, 7r/4) 


0.921405 
0.921707 


0.905433 
0.905979 


0.902971 
0.902974 


0.894779 
0.894999 


0.892836 
0.892828 


0.887437 
0.887552 


0.885867 
0.885852 


PNum(«, 7I"/2) 
PAsym(«, 7t/2) 


0.731659 
0.734432 


0.658904 
0.659884 


0.671538 
0.672036 


0.631168 
0.631562 


0.640067 
0.640252 


0.612185 
0.612384 


0.619215 
0.619276 


/0Num(™, 37r/4) 
PAsym(n,37r/4) 


0.541914 
0.548200 


0.338150 
0.338908 


0.444087 
0.445270 


0.312623 
0.312938 


0.395710 
0.396134 


0.295270 
0.295431 


0.365120 
0.365265 


Table A. 4: Numerical values of p(n, 9) for A = 0.3 


n 


2 


3 


4 


5 


6 


7 


8 


PNum(™,7r/4) 
PAsym(«, 7T/4) 


0.932056 
0.931684 


0.912068 
0.913012 


0.915471 
0.915038 


0.904519 
0.905020 


0.906536 
0.906242 


0.899088 
0.899399 


0.900492 
0.900275 


PNum(«, 7T/2) 
PAsym(«, 7T/2) 


0.768025 
0.768824 


0.676242 
0.677927 


0.712563 
0.712076 


0.656541 
0.657389 


0.683506 
0.683096 


0.642401 
0.642903 


0.664320 
0.663983 


/9Num("-,37r/4) 
PAsym(n,37T/4) 


0.603994 
0.607173 


0.354168 
0.355419 


0.511302 
0.511188 


0.335990 
0.336608 


0.464170 
0.463845 


0.322975 
0.323332 


0.433890 
0.433544 


Table A. 5: Numerical values of p(n, 9) for A = 0.7 


n 


2 


3 


4 


5 


6 


7 


8 


PNum(«, 7T/4) 
PAsym(n, 7T/4) 


0.903373 
0.903927 


0.887893 
0.887947 


0.878937 
0.879042 


0.870856 
0.870896 


0.865109 
0.865149 


0.859693 
0.859714 


0.855477 
0.855487 


/ONum(«,7r/2) 
PAsym(«, 7T/2) 


0.670096 
0.674020 


0.61307 
0.613330 


0.597811 
0.598582 


0.569629 
0.569822 


0.559974 
0.560266 


0.541915 
0.542024 


0.534903 
0.535013 


PNum(n,37T/4) 
PAsym(n,37T/4) 


0.436818 
0.446622 


0.295804 
0.296198 


0.332446 
0.334401 


0.256666 
0.256936 


0.283712 
0.284437 


0.232404 
0.232560 


0.253962 
0.254249 



Table A. 6: Numerical values of p(n, 9) for A = —0.3 
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n 


2 


3 


4 


5 


6 


7 


8 


/ONum(n,7r/4) 

PAsym(«, 7T/4) 


0.895942 
0.895477 


0.877985 
0.877611 


0.866417 
0.866290 


0.857132 
0.857039 


0.849944 
0.849893 


0.843718 
0.843674 


0.838522 
0.838481 


PNum(«, 7T/2) 
PAsym(«, 7T/2) 


0.644723 
0.646724 


0.587178 
0.586707 


0.562522 
0.562879 


0.535287 
0.535292 


0.520367 
0.520512 


0.503241 
0.503266 


0.492689 
0.492719 


/0Num("-,37r/4) 
/3A S ym("-,37T/4) 


0.393503 
0.402915 


0.271884 
0.271927 


0.284932 
0.286782 


0.226344 
0.226578 


0.236472 
0.237163 


0.199498 
0.199648 


0.207607 
0.207867 



Table A. 7: Numerical values of p(n, 9) for A = —0.5 



n 


2 


3 


4 


5 


6 


7 


8 


PNum(n,7r/4) 

PAsym(n, 7r/4) 


0.886810 
0.883755 


0.863491 
0.861533 


0.847596 
0.846496 


0.835540 
0.834864 


0.825990 
0.825552 


0.818029 
0.817718 


0.811263 
0.811018 


PNum(«, 7T/2) 
PAsym(n, 7T/2) 


0.613546 
0.610618 


0.549304 
0.545990 


0.512895 
0.511485 


0.483364 
0.482427 


0.462734 
0.462241 


0.444631 
0.444263 


0.430668 
0.430380 


/9Num("-,37T/4) 
PAsym("-,37T/4) 


0.340281 
0.348075 


0.236892 
0.235495 


0.225198 
0.226127 


0.182532 
0.182395 


0.177268 
0.177655 


0.153088 
0.153117 


0.150087 
0.150189 



Table A. 8: Numerical values of p(n, 9) for A = -0.7 



Appendix B Proof of (3.6) 

We prove Eq. (3.6) in the form 



4vr 2 

e^o 9 2 



2-2r? 



l/T) 



exp 



Jo 



sinh((27? - l)t) 2rj - 1 _ 2t \ dt 
sinhr/i cosh((l — rj)t) rj J t 



(B.l) 



By taking the logarithm and introduce a variable z = 0/2n, we can calculate the LHS from the 
definition Eq. (3.3) as follows. 



4tt 2 

limlnl D{2k-9)— r 
e^o V 9 2 



= In 



r 



(2-2*7) 



1/7, 



— lim 



sinh 2 (l - z)t 



(l-z)\_ 2t \ dt 



, , e " I • 2 In • 

z ^°[Jo V smn * c osh(l — 7])tsmhr]t 77 ) t 



Now substitute the function In z by its integral represention 

lnz= / (. 
Jo 



At r°° r]f 



e'-e -,-= / (e- 2t -e~^)^, 



(R(z) > 0) 



(B.2) 



(B.3) 
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we have 



2 



. , . sinh 2 (l - z)t (1 - z) 2 2 A dt 

lim { I [ — — ^ ^ ^ J -e~ 2t ) — + 2\nz 

sinnt cosh(l — rj)t smhrjt r] J t 



lim ( f ( . -^'"-f „ - 2e- + f 2 - <I^A .-»') * 
z ~*° [Jo \ sinnt cosh(l — rf)t sinhryt \ V J J t 

sinhi 2r? — 1 9 A dt 
+ 2 i e~ 2 * — 



/o \cosh(l — rj)t sinhryt 77 / t 

= f°° ( sinh((2r ? -l)t) _ 2r7-l c _ 2 A dt 
7 \ cos h(l ~~ r /)* sinhryt rj J t 



Thus Eq.(B.l), namely, Eq.(3.6) is proved. 



Appendix C Proof of (3.22) 



If we use the notation z = 0/(2w), the asymptotic amplitude of the string correlation function 
(3.3) is rewritten as 



2-2-q 



z 2 /ri 



cxp 



-/( 

J \smht 



sinh zt 



z 2 e~ 2t \ dt 



,cosh(l — r))t sinh rjt rj 



Setting the parameters as 

7 = it/2, A = cos7 = 0, rj = (tt — rf/ir = 1/2, 

we obtain 



D(9,tt/2) 



2:- 



cxp 



sinh 2 (zi) 



= 4 z exp 



-2 



sinhtcosh(t/2)sinh(t/2) 
( sinh 2 (zt) z 2 e -2t\ &~ 



2z 2 e~ 2t ) - 



V sinh t 



t 



(C.l) 



(C.2) 



(C.3) 



On the other hand, the Barnes G-function enjoys the following integral representation [48] 



roc -t / 2 f 2 \ 

m G (i + ,)=^ ^__ ¥ (i-^ + _-e-)d/ 



(l+7B)^+(Iog^)f, (»(*)>-!) 



(C.4) 
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where 7 is the Euler-Mascheroni constant. From this integral representation we have 

/•oo -t 

- In [G(l + z)G{\ -z)}= / — ^ (-2 - z 2 t 2 + 2 cosh(zt)) dt + (1 + lE )z 2 

Jo t(l — e ) 

-2 - z 2 t 2 + 2cosh(zt) dt 2 



4sinh^(t/2) t 
Z" 00 -1 - 2z 2 t 2 + cosh(2zt) dt „ , 2 

sinh 2 (zt) - z 2 t 2 dt 



f 

Jo 



+ (1+je)z 2 . 



sinh 2 1 

Now substituting the formula 



-2t 1°° 

= 1, 





_e 2 *-l 

into the standard integral representation of the Euler-Mascheroni constant, 

we can establish another integral formula 

/■°° f * 2 -2A dt 

Substituting Eq. (C.8) into Eq. (C.5) yields 

, r^/ 1 m /"°° /sinh 2 (zt) 2 _ 2t \ dt 

- In [G(l + *)G(1 - z)\ = ^ ^ s . nh V 2 ^ - z 2 e 2 < ' 

from which we conclude 

15(0, tt/2) = 4~ z2 [G(l + z)G(l - z)] 2 , 
which is equivalent to Eq. (3.22). 
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